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Recently PAMELA released their first results on the positron and antiproton ra- 
tios. Stimulated by the new data, we studied the cosmic ray propagation models and 
calculated the secondary positron and antiproton spectra. The low energy positron 
ratio can be consistent with data in the convection propagation model. Above 10 
GeV PAMELA data shows a clear excess on the positron ratio. However, the sec- 
ondary antiproton is roughly consistent with data. The positron excess may be a 
direct evidence of dark matter annihilation or decay. We compare the positron and 
anti-proton spectra with data by assuming dark matter annihilates or decays into 
different final states. The PAMELA data actually excludes quark pairs being the 
main final states, disfavors gauge boson final states. Only in the case of leptonic 
final states the positron and anti-proton spectra can be explained simultaneously. 
We also compare the decaying and annihilating dark matter scenarios to account for 
the PAMELA results and prefer to the decaying dark matter. Finally we consider a 
decaying neutralino dark matter model in the frame of supersymmetry with R-parity 
violation. The PAMELA data is well fitted with neutralino mass 600 ~ 2000 GeV 
and life time ~ 10^^ seconds. We also demonstrate that neutralino with mass around 
2TeV can fit PAMELA and ATIC data simuhaneously. 
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PACS numbers: 13.15.+g, 95.35. +d, 95.55.Vj, 98.62. Gq 



I. INTRODUCTION 

The existence of dark matter (DM) has been confirmed by many astronomical observa- 
tions, but the nature of DM is still an open question. Many kinds of particles in theories 
beyond the standard model (SM) are proposed as DM candidates The DM particles are 
usually stable due to the protection of some discrete symmetry. For example, the lightest 
supersymmetric particle (LSP) in the usual supersymmetric (SUSY) model is stale due to 
the conservation of R-parity, and is a well motivated candidate for DM 2|. 

Generally there are three kinds of strategies to detect the DM particles: the "collider 
search" at large colliders such as large hadron collider (LHC) or international linear collider 
(ILC); the "direct detection" to find the signal of nuclei recoil when DM particles scatter 
off the detector; the "indirect detection" to search for the products from DM annihilation 
or decay, such as neutrinos, photons and anti-matter particles. It is usually a big challenge 
for indirect detection due to the difficulty in discriminating the signal from the astrophys- 
ical background. Therefore precise measurements with wide energy range and improved 
resolution are necessary for DM indirect detection. 

PAMELA is a satellite borne experiment designed to measure the cosmic rays (CRs) in a 
wide energy range with unprecedented accuracy Recently, the PAMELA collaboration 
released the first data about antiprotons and positrons % ISj]. Usually it is thought that 
antiprotons and positrons are produced when CRs propagate in the Milky Way and collide 
with the interstellar medium (ISM). The abundance of these secondaries can be calculated 
with relatively high precision. However, the PAMELA results show an obvious excess in the 
faction of e"^/(e"'" + e~) at energies above ~ 10 GeV. Interestingly the excess keeps to rise 
up to energy ~ 100 GeV. On the other hand, the spectrum of antiprotons fits the prediction 
quite well. These results confirm the previous results by HEAT jo] and AMS within the 
error bars. 

The PAMELA results may provide a direct evidence of DM in the way of "indirect detec- 
tion". However, before resorting to the exotic physics of DM, it is necessary to go through 
the possible astrophysical sources to account for this results. The model-independent spec- 
tral shape analysis shows that there might be most likely a primary source with e^e~ pairs 
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required to explain the rise of the positron fraction 



The non-excess of the antiproton data 



also favors a leptonic origin of the positrons. Pulsar is thought to be a good candidate to 
produce only leptons, and was used to explain the previous HEAT data [l^ • The recent 
analysis show that the PAMELA data can also be fitted considering the contributions from 
nearby pulsars such as Geminga and B0656+14 
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interesting conclusion in Ref. 



12|. It should be noted that another 



131 ] shows that the uncertainties of the propagation of CRs, 



the production cross section of secondary particles and the errors of electron measurements 
might lead to the underestimation of the positron fraction, and the "excess" is actually not 
an excess. 

As one possibility, the contribution from DM annihilation or decay is also widely dis- 
cussed. The scenarios about DM annihilation include (i) annihilation to SM particle pairs. 



like gauge boson, quark and lepton pairs 



process to e+e 
to e~^e 
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171 ] ; (ii) virtual internal bremsstrahlung 



; (iii) annihilation to new mediating particle pair which would decay 
etc. However, in the DM annihilation scenario, an unnatural large 
"boost factor" is necessary to reproduce the PAMELA positron data. Another problem is 
that the non-excess of the p/p data will set constraint on the properties of annihilating DM 



which makes the model building difficult. 
Considering the difficulties to explain the large positron excess by annihilating DM, we 
propose to solve the problem using decaying DM in the present work. If there exists some 
tiny symmetry violation, the DM particles can decay very slowly to SM particles. The 
decaying DM models are studied extensively in literatures to explain the observational data 
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34]. 
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291 ] or set constraints on the decay properties of DM [30 , 
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For the self-consistence of this work, we first studied the process of cosmic ray propagation 
carefully, and give realistic propagation models to produce the background contributions to 
positrons and antiprotons from CR interaction with ISM. The DM induced positrons and 
antiprotons are calculated in the same propagation models. We find that the PAMELA data 
can be well fitted for a leptonically decaying DM (LDDM) model. The results for different 
decay final states are discussed in detail. And a possible model in the SUSY frame with 
R-parity violation is proposed. 

The paper is organized as follows. In section II, we discuss the propagation of CRs and 
give the updated positron and antiproton background estimations. In section III, we present 
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a model-independent approach to recover the PAMELA data and discuss why we need the 
decaying DM whose decay products are mainly leptons. In section IV, we give an example 
of this LDDM model and discuss another possibility. Finally we give the summary and 
discussion in section V. 



II. PROPAGATION OF GALACTIC COSMIC RAYS 



In this section we will study the cosmic ray propagation model carefully so that we 
can predict the positron and antiproton spectra to compare with the PAMELA data. The 
charged particles propagate diffusively in the Galaxy due to the scattering with random 



magnetic field 



35|. The interactions with ISM and interstellar radiation field (ISRF) will 



lead to energy losses of CRs. For heavy nuclei and unstable nuclei there are fragmentation 
processes by collisions with ISM and radiactive decays respectively. In addition, the overall 
convection driven by the Galactic wind and reacceleration due to the interstellar shock will 



also affect the distribution function of CRs. The propagation equation can be written as 



^ = Q(x,p) + V . (B,,VV' - V,^) + -■p'D„--^,p 

-^U-f(v.vw-)l-^-^. (1) 

op I 3 \ Tf Tr 

where if) is the density of cosmic ray particles per unit momentum interval, (5(x,p) is the 
source term, D^x is the spatial diffusion coefficient, Vc is the convection velocity, Dpp is 
the diffusion coefficient in momentum space used to describe the reacceleration process, 
p = dp/dt is the momentum loss rate, rj and are time scales for fragmentation and 
radioactive decay respectively. We describe a bit more about the relevant terms in Eq.([T]) 
in the follow. 

For primary particles such as the protons and some heavy nuclei, the source function is 
the product of two parts: the spatial distribution /(x) and energy spectrum q{p). /(x) can 
bllow the distribution of possible sources of CRs, such as the supernova remnants (SNR) 
371]. The injection spectrum oc p~'^ is usually assumed to be a power law or broken 
power law function with respect to momentum p. For secondary particles the source function 
is given according to the distributions of primary CRs and ISM 



Q(x,p) = /3c?/^p(x,p)[crH(p)nH(x) + CrHc(p)?^Hc(x)], 



(2) 
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where il)p{'x,p) is the density of primary CRs, f3c is the velocity of injection CRs, (Th and ane 
are the cross sections for the secondary particles from the progenitors of H and He targets, 
riH and riHe are the interstellar Hydrogen and Helium number densities, respectively. 

The spatial diffusion is regarded as isotropic and described using a rigidity dependent 
function 

D.. = f3Do 

The reacceleration is described by the diffusion in momentum space. The momentum diffu 



(3) 



sion coefficient Dpp relates with the spatial diffusion coefficient Dr^x as 



38| 



DppDxx 



(4) 



35(4-(52)(4-5)w' 

where va is the Alfven speed, w is the ratio of magnetohydrodynamic wave energy density 
to the magnetic field energy density, which characterizes the level of turbulence, w can be 



taken as 1 and the reacceleration is determined by the Alfven speed va i 

The convection velocity, which corresponds to the Galactic wind, is assumed to be cylin- 
drically symmetric and increase linearly with the height z from the Galactic plane 36|] . It 
means a constant adiabatic energy loss of CRs. Vc{z = 0) = is adopted to avoid the 
discontinuity across the Galactic plane. 

Finally the energy losses and fragmentations can be calculated according to the interac- 
tions between CRs and ISM or ISRF. 



For some simplified cases t 
the Green's function method 



le p rop agation equation ([T]) can be solved analytically using 
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4l|. However, generally it is not easy to find the 



analytical solution. A numerical method to solve this equation has been developed by Strong 



and Moskalenko, known as the GALPROP model [36|, 142]. In GALPROP, the realistic 
astrophysical inputs such as the ISM and ISRF are adopted to calculate the fragmentations 
and energy losses of CRs. The parameters are tuned to reproduce the observational CR 
spectra at Earth. It is shown that the GALPROP model can give relative good descriptions 



of all kinds of CRs, including the secondaries such as e"*", p and diffuse 7-rays 
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45|. 



In this work we employ GALPROP models to calculate the propagation of CRs. Two 
GALPROP models are adopted. One is the diffusion + convection (DC) scenario and the 



other is the diffusion + reacceleration (DR) model. It has been shown t 
is easier to reproduce the energy dependence of the observed B/C data[i 



lat the DR model 



36|. However, the 
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TABLE I: The propagation parameters in the DC and DR models. 



^0 


diffusion index" va dVc/dz 


e injection* 


nuclei injection'^ 


(10^8 cm2 s-^) 


81/62 (km s~"^) (km kpc~^) 


71/72 


71/72 


DC 2.5 


0/0.55 — 6 


1.50/2.54 


2.45/2.25 


DR 5.5 


0.34/0.34 32 — 


1.50/2.54 


1.94/2.42 



"Below/above the break rigidity po = 4 GV. 
''Below/above 4 GcV. 

'^The break energy is 25 GeV for DC model, and 15 GeV for DR model. 



reacceleration will produce more low energy CRs and overestimate the low energy spectra 



of electrons, positrons, protons and Helium 46j. In addition the DR model seems to un- 
derproduce the antiprotons 43|, |47l]. In DC model, the results of e~, e"*", p and p are in 
better agreement with the data, but the "peak" around 1 GeV of the B/C data is not well 
generated in the model 36|. Here we quote these two models in the sense that the differences 
between these models are regarded as the uncertainties of the propagation model of CRs. 
The propagation parameters are listed in Table [B Other parameters which are not included 
in the table are: the height of the propagation halo 2/^ = 4 kpc, the spatial distribution of 



primary CRs /(x) oc 



R 

^0 



exp 



/3(R-Rq) 
Rr 



exp 



with a = 0.5, p = 1.0, R, 



8.5 



kpc and Zg = 0.2 kpc 481]. 

In FigHlwe show the observed and calculated CR spectra of B/C, ^'^Be/^Be, protons and 
electrons for both the DC and DR models. For the solar modulation of the local interstellar 
(LIS) spectrum we adopt the force field approximation [gS]. It is shown that these two models 
can both give satisfactory descriptions of the data. We can also note that the DR model 
indeed produce more low energy electrons, and a larger solar modulation potential is needed 
to suppress the low energy spectra. 

Fig. [2] gives the results of positron fraction and p/p ratio for DC and DR model re- 
spectively. We can see from this figure that the DR model gives too many positrons at 
low energies. The solar modulation does not change the results significantly since it affects 
positrons and electrons simultaneously. The charge dependent solar modulation effect might 
be helpful in softening the discrepancy between the calculation and data. The result of DC 
model shows a better agreement with the low energy data. We will not focus on the low 
energy behavior of the positron fraction since it might be mainly due to the solar effect. For 
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FIG. 1: Propagated B/C, ^"Be/^Be, protons and electrons spectra in the DC and DR mod- 
els of GALPROP. In each panel of the figure, the thick solid lines represent the results of DC 
model, while the thin solid lines show results of DR model. For each model the LIS spectrum 
together with the solar modulated one are plotted. In order to match the low energy data, dif- 
ferent modulation potentials are adop ted as labeled in the figure. The references of the data are, 
B/C: Chapell-Web 
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59|; protons: 



energies higher than several GeV, the results of the two models are similar, and both seem 
to underestimate the positron fraction compared with the HEAT jol and PAMELA data . 
As for p/p, DC model is in good consistent with the measurements, including the recent 
PAMELA dataj4|. The DR model shows an underproduction of antiprotons. It means that 
if the DC model is correct, the excess of positrons and non-excess of antiprotons will set 
strong constraints on the properties of the source of positrons, e.g., l^j- On the other hand, 
the DR model will leave more loose constraint. 
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FIG. 2: Left: the calculated positron fraction compared with observations; right: p/p ra- 



tio. References of the observational data are, positron fraction: TS93 
fl, HEAT94+95 3], HEATOO 67], PAMELA |3; p/p: IMAX P 
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III. MOTIVATION FOR A LEPTONICALLY DECAYING DM MODEL 

In the section we will adopt a model independent approach to constrain the DM anni- 
hilation or decay products from the PAMELA data. We find the PAMELA data actually 
excludes the annihilation or decay products being quark pairs, strongly disfavors the gauge 
bosons and favors dominant leptonic final states. 

In principle, there should be no difference in treating annihilation or decay. However, 
there is a subtle difference for the two scenarios. We know that the annihilation or decay 
depends on the DM density in different ways. Further since what we observed on Earth 
is the integrated positron or antiproton flux from the nearby region, we may get slightly 
different spectra after propagation for the annihilation or decay scenarios even the source 
spectra are the same. In the section we only study the case of decay to show the result. 
However, the main conclusion will be unchanged for annihilation. In the last of the section 
we also show why decaying DM is superior to annihilating DM. 

The source term in Eq. [1] is given by 

V ~ \decay \p) 

TDM rUDM dE 

where p(r) is the DM density distribution in the Galaxy, tdm is fhe life time of DM, dN/dE 
is the original positrons spectrum from each DM decay. The DM mass niDM is supposed to 



9 



be a free parameter while its life t^m will be fixed by the positron fraction. In the work we 
take the NFW profile for DM distribution with the local DM density Pq = 0.3GeV/cm^. 

We assume two-body final state with energy of 100 and 300 GeV for either quark, lepton 
or gauge boson pairs (1 TeV is also discussed in this case). The spectra of positron and 
antiproton are simulated by PYTHIA 75|] and then propagated by adopting the DC and 
DR propagation models as introduced in the last section. The life time of DM is taken to 
give correct positron flux. 




-bg 

-w100GeV1.5 
■z 100GeV1.1 
-w300GeV1.6 
-z 300GeV1.5 
-w1000GeV1.6 
-z 1000GeV1.l 



10 




bg 

w100GeV1.4 

z 100GCV1.2 

w300GeV1.2 

- - z 300GeV1.2 

w 1000GeV1.2 

z 1000GeV1.2 



(GeV) 



1E-3 



a 
ia 



1E-S 



1 

DC 


M 1 ' 
















r# — >^ 






-— "^^^ X jfr 










bg 




w100GeV1.5 




z 100GeV1.1 - 




w300GeV1.6 : 




- - z 300GeV1.5 '■ 




w1000GeV1.6 ' 


' 


z 1000GeV1.6 ' 





0.1 



10 

, (GeV) 



100 



1E-3 



a 
ia 



1E-S 



1 1 

- DR 


1 ' 
































bg 




w100GeV1.4 




z 100GeV1.2 ■; 




w300GeV1.2 i 




z 300GeV1.2 ' 




w1000GeV1.2 ' 





z 1000GeV1.2 ' 





0.1 



10 

, (GeV) 



100 



FIG. 3: The positron and antiproton fraction, ^^^^ and as function of energy from DM 
decaying to gauge boson pairs. The black lines are background for positron and anti-proton. The 
data points are from the preliminary PAMELA results J]. The numbers 100, 300 and 1000 GeV 
refer to the energy of gauge bosons, while the numbers 3.6 and so on refer to the life time of DM 
in unit of 10^^ s. In the title, "DC" and "DR" are the two propagation models discussed in the 
section II. 
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In FigJ3l we show the positron and antiproton fraction when the decay products are gauge 
boson pairs. We find that the positron spectrum from gauge boson decay is usually softer 
than PAMELA given even taking the gauge boson energy at 1000 GeV. Especially, we find 
the gauge boson channel is problematic for the antiproton spectrum. They give several times 
larger anti-proton fraction than the data in the two transportation models. Therefore DM 



decaying to gauge 



Dosons are strongly disfavored by the anti-proton data. This conclusion 



is similar with Ref . l3, H] the authors suggest to use extremely high energy gauge bosons of 
about lOTeV to interpret both positron and antiproton spectra. In that case, only the soft 
tail of positron and antiproton spectra from high energy gauge boson are adopted. Since the 
endpoint of the anti-proton soft tail has energy larger than 100 GeV it has no confiict with 
the present PAMELA data that cut off at ~ 100 GeV. It predicts very high antiproton flux 
above ~ 100 GeV. The model requires very heavy DM (~ 10 TeV). It should be mentioned 



that the PPB-BETS 



and ATIClZSl] see electron spectrum has a 

signiflcant excess above several hundred GeV with a cutoff around SOOGeV. The heavy DM 
with mass of lOTeV seems not to account for this positron excess 17 1. 

In FiglU we show the positron and antiproton fraction when DM decay to quark pairs. 
Positrons are produced after hadronization of quarks via the decay of charged pions. However 
we flnd these positrons are too soft and can not account for the excess of positrons above 
~ 10 GeV. Furthermore, quark hadronization produces unwanted anti-protons as expected 
which are several times larger than the experimental data. 

In Figini we show the case of leptonic flnal states. The positron spectrum is easily to 
flt the PAMELA data in this case for e and /i flnal states. However, for r flnal states we 
may need the r energy larger than ~ 300 GeV to account for the hard positron spectrum. 
Certainly there is negligible influence on the antiproton spectrum in the pure leptonic decay. 





Gauge boson 


Quarks 


Leptons 


Positron 




X 




anti-Proton 


X 


X 





TABLE II: The summary of three kinds of decay products to account for the PAMELA data. 

In TableJIT] we give a summary on these three channels to flt the PAMELA data. The 
gauge bosons flnal states with energy around hundreds of GeV to ~TeV have problems 
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FIG. 4: Positron and antiproton fraction as function of energy from DM decaying to quark pairs. 
The labels are the same as that in Fig. [3l 

with anti-protons ^. However if the energy is extremely high (~ 10 TeV), the positron and 
anti-proton data can be satisfied due to the energy cutoff of the PAMELA data. The quark 
final states have problem in both the positron and anti-proton spectra, that is they give too 
soft positron spectrum and to large antiproton flux. The lepton final states with hundred 
GeV can give hard positron spectrum and easily fit the PAMELA data very well. At the 
same time they are free from upsetting the anti-proton spectrum. 

Finally we will try to compare the two scenarios of annihilating and decaying DM. We 
find, if the positron excess at PAMELA is indeed of DM origin, the decaying DM is superior 
to annihilating DM. First, annihilating DM has to resort to large "boost factor" at the order 

^ The conclusion is drawn in the conventional propagation model. It is possible to moderate or overcome 
the problem in some special propagation models. 
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FIG. 5: Positron and antiproton fraction as function of energy from DM decaying to lepton pairs. 
The labels are the same as that in Fig. [3l 

of ~ 10^ to ~ 10^ to account for the large positron flux. However, detailed analysis of boost 
factor from the dumpiness of DM structures based on N-body simulation gives that the 
most probable boost factor should be less than 10 — 20 3l- The same conclusion is also 



found through the direct computation of the antimatter fluxes from N-body simulation 801]. 
A nearby subhalo or the DM spike around the intermediate mass black hole might be able to 



provide large boost factors, however, these scenarios are found to be of little probability 
or suffer large uncertainties 



821. Some authors use Sommerfeld-enhancement to increase 



the annihilation cross-section 



ill 



14 



17 



19 



20 



83 



84| and some people adopt non-thermal 



DM such as wino 85|].^ This method is usually constrained by the data of gamma rays 
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881] and need further study. Second, annihilating DM usually produce more quarks or 



gauge bosons than leptonic final states. For example neutralino is more easily to annihilate 
into quarks or gauge bosons than leptons. The Kaluza-Klein DM in some universal extra 
dimension models may be an exception that it can annihilate largely to leptons. However, we 
should mention that Kaluza-Klein DM may still have problem on the PAMELA antiproton 
data, because in its annihilation final states quarks and gauge bosons are comparable to 



^ It is interesting to mention that DM may be more than one-component, for example, one component is 
metastablc in the early Universe and the other component is stable that it can annihilate into leptons 
today which can also avoid large boost factor [s" 
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leptons. From Figs. El IH and we can roughly estimate that they may still produce too 
much antiprotons. 

In the next section, we will consider a scenario of LDDM. Our example is given in the 
minimal supersymmetric standard model with trilinear R-parity violation term LLE. The 
neutralino being the lightest SUSY particle forms dark matter and is produced thermally at 
early universe in the same way as stable neutralino. To account for the PAMELA positron 
excess we find it has the life time of about ~ 10^^ seconds, which is much larger than the 
life of the universe. 



IV. NEUTRALINO WITH R-PARITY VIOLATION AS AN EXAMPLE OF 

LDDM 



In this section we turn to discuss the specific decaying DM model in supersymmetry 
scenario. In SUSY model, the discrete symmetry of R-parity can be invoked to avoid dan- 
gerous baryon-number violation terms which drive unexpected proton decay. Defined as 
R = (_i)2S'+3B+L^ ^Yie R-parity of a SM particle is even while its superpartner is odd. 
Then the LSP particle is stable. Since the neutralino LSP can have correct relic density via 
thermal production which makes it a suitable DM candidate. 

Although R-parity symmetry is well motivated for SUSY phenomenology, there is no 
reason for this symmetry to be exact. One can introduce some R-parity violation terms in 
the Lagrangian which make the LSP decay into SM particles. The general gauge invariant 
superpotential of the minimal supersymmetric standard model can be written as 



W = Wmssm + XijkLiLjEk + X'ij^LiQjDk + X'-j^UiDjDk + jiiLiHu (6) 

where i, j, k are generation indices (we neglect these indices bellows). The LH term in this 
superpotential mixes the lepton and Higgs field. Because Higgs and Lepton have the same 
gauge quantum numbers, one can rotate away the bilinear term by redefining these fields. 
The LLE, LQD terms violate lepton number, and the UDD terms violate baryon number 



89[. Then totally 45 couplings including 9 Xijk, STA^^;. and OA-'^-^, are invoked in theory. 
There might also exist many soft SUSY breaking terms which induce R-parity symmetry 
breaking. These terms added more additional free parameters, and we do not discuss them 
here. 
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The R-parity violation terms must be tiny to satisfy stringent experiment constraints, 
especially the constraints from proton decay. The tiny R-parity violation may be due to 
some fundamental theories at high energy scale. For example, in SU(5) grand unified the- 
ories (GUT) the gauge invariant terms as fijk5^5^10^ , where 5, 10 denote matter field rep- 



resentations contain leptons and quarks, could induce R 



90, 



)ariw 
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violation explicitly, fijk is 



92l ] . To construct a realizable 



suppressed and gives tiny R-parity violation couplings 
theory taking account for all the experimental constraints, more symmetries and high order 
operators are always required. 

In general, we should consider all R-parity violation terms in Eq. |6l However, they might 
not appear in theory at the same time. In some theories only baryon-number violation terms 
UDD is include d | 93l| . while some other theories can only predict lepton- number violation 



LQD term 
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In a class of discrete gauge symmetric models the baryon-number 



9J] . In this work. 



violation term should be absent since they may induce rapid proton decay 
we only consider lepton-number violation term LLE in order to account for the PAMELA 
positron excess. We assume all the components of are equal and neglect the annihilation 
of neutralinos in the following discussion for simplicity. 





SUSY 


MC 


Mass(GeV) 


mo(GeU) 


mii2{GeV) 




tan P 


A 


SPS6 


bino 


190 


150 


300 


1.04 X 10"^ 


10 




SUSY 


MC 


Mass(GeV) 


mo(GeU) 


mi/2(GeU) 


Qh^ 


tan P 


B 


mSUGRA 


bino 


341 


900 


800 


9.62 X 10"2 


50 


C 


mSUGRA 


bino 


614 


1750 


1400 


9.97 X 10-2 


50 


D 


mSUGRA 


bino 


899 


5000 


2000 


1.02 X 10"^ 


50 


E 


mSUGRA 


higgsino 


1126 


9100 


3500 


1.01 X 10"^ 


50 




SUSY 


MC 


Mass(GeV) 


mo(GeU) 


m3/2(GeU) 




tan P 


F 


AMSB 


wino 


2040 


18000 


640000 


9.15 X 10^2 


10 


G 


AMSB 


wino 


2319 


20000 


730000 


1.17 X 10"^ 


10 



TABLE III: The seven benchmark points with different neutralino masses and in different scenarios. 
"MC" means the main component of neutralino. The other parameters we adopted are sgn(/i) = 



+l,mt = 172.6GeV,Ao = 0, 
Ms at GUT scale. Suspect 



he SPS6 point has non-universal mass parameter Mi/1.6 = M2 



95l | and MicrOmega 



96l ] are used in our calculation. 
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To specify the SUSY parameters we choose some benchmark points as denoted in Tab. 
Illlt where LSP is neutrahno with different masses. The point A is SPS 6 in the mSUGRA- 
hke scenario with non-unified gaugino masses, the point B, C, D are from the Higgs funnel 
region in the mSUGRA, the point E is from focus point region in the mSUGRA, and the 
point F denotes a thermal wino in the AMSB scenario. All these points satisfy the correct 
relic density and other laboratory constraints. We utilize PYTHIA to produce the positron 
energy spectrum from neutralino decay. 

In the limit of heavy and degenerate sleptons, neutralino with gaugino component has 



life time 



32| 



V 10-25 y yiOOOGeVJ 

icient for LLt 
higgsino component has life time 



where A' is the coefficient for LLe term and is the sfermion mass. The neutralino with 



^26„ /tan/^yV ^' VV Y'f^r" 



DC 


r(1026s) 


A'(10-25) 


DR 


t(1026s) 


A' (10-25) 


A 


4.5 


3.2 


A 


3.9 


3.4 


B 


2.6 


14.8 


B 


2.3 


15.7 


C 


1.7 


16.1 


C 


1.5 


17.1 


D 


1.2 


59.5 


D 


1.1 


62.1 


E 


1.0 


253.2 


E 


0.9 


266.9 


F 


0.6 


159.7 


F 


0.6 


159.7 


G 


0.5 


156.7 


G 


0.5 


156.7 



TABLE IV: Life time r in unit of lO^^s and R-parity violation parameter A' in unit 10 ^5 for 
different benchmark points. "DC" and "DR" refer to the two propagation models we adopted. 

These neutralinos generally have life time around lO^^s to account for the PAMELA data. 
It is much longer than the life of the universe which makes these neutralinos as valid dark 
matter. However, the tiny value of A' means we can not test the scenario directly at LHC 
or ILC. In Tab. [IV]we give the life time and value of A' for the different benchmark points 
to account for the PAMELA positron excess. 
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E«„ (GeV) E^,„ (GeV) 



FIG. 6: Positron and antiproton fraction as function of energy for different benchmark points. 
The neutrahno mass of each model is given, while the life time of neutralino in each model is given 
in Tab llVl The neutralino mass in each line is increasing from left to right. 

In Fig. El we show the positron fraction including contribution from neutralino decay for 
different benchmark points. From Fig. [6l we can see the positrons are most probably dis- 
tributed lower than ^m^o. This is due to the three-body decay and the roughly equal decay 
fraction of Viljl^ ■ The positron spectrum given here with three-body decay is generally 
softer than two body final states with monochromatic lepton. We notice that neutralino 
heavier than 300GeV to about 2 TeV can fit the PAMELA data well. If one takes account 
the results from PPB-BETS and ATIC experiments, it seems to suggest a heavy neutralinos 
up to ~TeV if they are due to DM decay. Another feature of this scenario is neutrino flux 
from neutralino decay, but we find this neutrino flux is too small to be detected by neutrino 
telescope such as Super-kamiokande, IceCube, etc while which is well within the atmospheric 
neutrino bound from Super-kamiokande(similar analysis for neutrino flux from decaying DM 



can be found in Ref. 



27 



m 



Finally we will give some comments on other scenarios of decaying dark matter. In some 
SUSY scenarios, the LSP may be the neutral particles as gravitino or sneutrino. They 
can also be the candidates of DM. However, in the bilinear R-parity violation scenario, the 



gravitino would mainly decay into gauge boson plus lepton 



25 



26| which is not favored 



by antiproton data from PAMELA. If the R-parity violation term is LLE as we discussed 
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FIG. 7: Positron and electron total flux as function of energy for different benchmark points 
comparing with ATIC data Ref. 7g|. The labels are the same as that in Fig. [6l 



above, the spectrum of decay products is similar. But the decay of gravitino is suppressed 
by the Planck scale and the R-parity violation parameter can be as large as A^^^ > 10~^^, 



which may lead to unstable NLSP decay quickly to avoid 



destroying the success of big-bang 



nucleosynthesis 22, ISSj and be observable at colliders 22 1 . 



Another possible candidate of DM in SUSY is sneutrino. The left-handed sneutrino has 
been ruled out by DM direct detection, therefore the right-handed sneutrino receives more 
concerns. If we neglect all the soft breaking terms, the interactions for right-handed neutrino 
superfeild N are induced only from Yukawa term yijLiHuNj. After introducing only bilinear 
R-parity violation terms ^[LiHu, the right-handed sneutrino will decay into lepton pairs with 
a large fraction, due to mixing between Higgs and left-handed lepton superfields. However, 
it should be mentioned that mixing between right-handed sneutrino and Higgs boson can 
not be avoided in general. These mixing terms have to be suppressed in order not to conflict 
with PAMELA antiproton data. The kinematical conditions can be used to suppress the 
right-handed sneutrino decay to two Higgs or one Higgs plus one lepton. To suppress the 
LSP direct decay to quarks which is induced by the mixing of right-handed sneutrino and 
Higgs, an extra relation Uiju'^ ~ is needed for the j-th generation right-handed sneutrino 



which is the LSP. In the Ref. 



29|] authors proposed similar idea, they assume the LSP is the 



third generation right-handed sneutrino and only /i'-^ ^ 0, then the small yukawa coupling 
yi3 make this right-handed sneutrino as a good candidate of LDDM. The difference between 
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leptonically decaying neutralino and right-handed sneutrino might be represented in the 
original positron spectrum. The former one is softer than the latter, because the former one 
is three-body decay while the latter one is two-body decay. It means that the PAMELA 
implication on DM mass is different for these two scenarios. With the synergy of PAMELA 
and collider experiments, it is possible to distinguish whether the DM is neutralino, gravitino 
or right-handed sneutrino. 

V. SUMMARY AND DISCUSSION 

The recent released PAMELA data shows interesting features, that the positron shows 
an obvious excess above ~ 10 GeV while the antiproton flux is well consistent with the 
expectation from the conventional cosmic ray model. The result implies there should exist 
some kinds of primary positron sources in addition to the secondaries from CRs interactions 
with the ISM. 

Before the discussion of possible exotic sources of positrons it is extremely important to 
explore the background carefully. In this work we first recalculate the background contribu- 
tions of positron and antiprotons from CRs using GALPROP. Considering the CR data of 
unstable secondaries, such as ^'^Be/^Be, it is possible to reduce the uncertainties in predicting 
the positron and antiproton flux from the propagation parameters |l^. Two typical propa- 
gation models, the diffusion-|-convection and diffusion-|-reacceleration ones, are adopted as 
benchmarks of the CR propagation models. In both of the two models, propagation param- 
eters are adjusted so that most of the CR spectra, such as B/C, ^°Be/^Be and so on, are 
well consistent with the observation. The DC model is found to be consistent with the p/p 
and the low energy /{e^ + data of PAMELA without introducing a charge- dependent 
solar modulation model. While the DR model shows a bit of underestimation of the p/p 
data and produce more low energy positrons. In both models the positron fraction at high 
energies shows obvious excess. 

We then consider DM origin of the primary positrons. We compared the positron and 
antiproton spectra with PAMELA data by assuming the annihilation or decay products 
are mainly gauge bosons, quarks and leptons respectively. We find the PAMELA data 
exclude quarks being dominant final states, disfavor gauge bosons and favor the leptonic 
final states. Comparing annihilating DM and decaying DM scenarios, we prefer decaying DM 
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because annihilating DM usually requires very large boost factors. Moreover, annihilating 
DM usually produce more gauge bosons and quarks than leptons. 

A concrete example of such LDDM model is considered in MSSM with tiny LLE R- 
parity violation term. We choose neutralino as the LSP and the DM particle. We show 
that neutralino with mass 600 ~ 2000 GeV and life time of ~ 10^^ seconds can fit the 
PAMELA data very well. We also demonstrate that neutralino with mass around 2TeV 
can fit PAMELA and ATIC data simultaneously. Another LDDM model is right-handed 
sneutrino being DM with bilinear R-parity violation term. With the interplay of PAMELA 
and LHC it is possible to distinguish the DM between neutralino and right-handed sneutrino. 

Finally we will discuss the implications of the LDDM models. LDDM can have many 
interesting phenomena in experiments other than PAMELA, such as gamma ray observation. 
In general, it can have two-body decay and three-body decay channels as DM — s> and 
DM I'^l'^X where the leptons in the final-states are required from PAMELA data. 
Particle X can be photons 7, neutrino v in general. If the LDDM has the decay channel 
DM — > l~^i~'~f, the spectrum of gamma ray can reveal the property of LDDM and the 
mechanics behind. For example, if there is an intermediating resonance particle Y which 
decays to lepton pairs, the decay of LDDM is now DM '-/Y l^l~'~f. Thus the photon is 
monochromatic and can be explored by the Fermi Gamma Ray Space Telescope (FGRST). 
Detecting the gamma ray spectrum will be an important way to test the LDDM model, 
where FGRST will play an important role. 

For our neutralino LDDM model the neutrino observation will provide a way to test the 
model since the decay produces hard neutrino spectrum. Another possible way to test the 
model is to look for the synchrotron emission of the hard electrons and positrons by DM 
decay. Further studies on this issue are on going. 
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